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Since the 1950s Heisenberg and others have attempted to explain the appear- 
ance of countable particles in quantum field theory in terms of stable localized 
field configurations As an exception Skyrme's model succeeded to describe 
nuclear particles as localized states, so-called 'skyrmions', within a non-linear 
field theory j^. Skyrmions are a characteristic of non-linear continuum models 
ranging from microscopic to cosmological scales j^, Q, ^, Skyrmionic states 
have been found under non-equilibrium conditions, or when stabilised by ex- 
ternal fields or the proliferation of topological defects. Examples are Turing 
patterns in classical liquids j^, spin textures in quantum Hall magnets (sl, or 
the blue phases in liquid crystals [9], respectively. However, it is believed that 
skyrmions cannot form spontaneous ground states like ferromagnetic or antifer- 
romagnetic order in magnetic materials. Here, we show theoretically that this 
assumption is wrong and that skyrmion textures may form spontaneously in con- 
densed matter systems with chiral interactions without the assistance of external 
fields or the proliferation of defects. We show this within a phenomenological 
continuum model, that is based on a few material-specific parameters that may 
be determined from experiment. As a new condition not considered before, we 
allow for softened amplitude variations of the magnetisation - a key property of, 
for instance, metallic magnets. Our model implies that spontaneous skyrmion 
lattice ground states may exist quite generally in a large number of materials, 
notably at surfaces and in thin films as well as in bulk compounds, where a lack 
of space inversion symmetry leads to chiral interactions. 

The possibility that particle-like states may form spontaneously in continuous fields has 
motivated intense theoretical efforts in the past. Derrick and Hobart established by rather 
general arguments that particle-like configurations are not stable in the majority of non- 
linear field models jlfll . However, a few exceptions have been found. Skyrme showed 



that particle-like excitations of continuous fields exist in the presence of certain non-linear 

n 

interactions 2]. As a drawback, the interactions considered by Skyrme are physically not 
transparent, because they involve higher order derivative terms that are technically in- 
tractable. Therefore, Skyrme's approach is not viable in the context of ordered states in 
condensed matter that are ruled by short range interactions. In contrast, a physically trans- 
parent exception to the Derrick-Hobart theorem has been recognized in systems with bro- 
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ken inversion symme try, where chiral interactions lead to skyrmion excitations in condensed 



matter systems 
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14 . llG^ . Chiral interactions exist in many different systems, e.g., (i) 
spin-orbit interactions in non-centrosymmetric materials, also referred to as Dzyaloshinsky- 
Moriya (DM) interactions 13], (ii) in non-centrosymmetric ferroelectrics, (iii) for certain 
structural phase transitions, (iv) in chiral liquid crystals, and (v) in the form of Chern- 
Simons terms in gauge field theories (see section V in |2Ci] for references). Yet, it was 
concluded that spontaneous skyrmion ground states are not stable, but can only be induced 
by an appropriate applied field j3| or the proliferation of topological defects Finally, a 
somewhat different mechanism has been found in quantum Hall ferromagnets. Here charge 
as an additional degree of freedom allows to stabilises skyrmions in an external field. As 
for the other examples skyrmions in quantum Hall ferromagnets do not form spontaneous 
ground states. 

We have chosen to study skyrmions that are driven by chiral interactions, because the 
magnetism of condensed matter systems provides perhaps the richest setting of physics 
problems in which different kinds of chiral interactions may be studied. Present day ex- 
perimental and theoretical studies have led to the view that only one-dimensional 
(uniform) modulations exist as spontaneous ground states in these systems (for an illus- 
tration see Fig. la). In the following, we show that previous studies have been incomplete 
and that skyrmions may form spontaneously in magnets with chiral interactions in terms of 
multi-dimensional modulations, when the amplitude of the magnetisation is exceptionally 
soft. The mathematical structure and the physical nature of the skyrmion texture consti- 
tute the first example for a spontaneous skyrmion ground state in an ordered condensed 
matter system. The example we identify is in particular important, because it can be tested 
experimentally in comparatively simple laboratory experiments. 

In our calculation we consider systems with a hierachy of well-separated energy and 
length scales, where ferromagnetic exchange favouring spin alignment is the strongest scale, 
followed by chiral interactions favouring spin rotations on intermediate scales, and mag- 
netic anisotropics and dipolar interactions that determine favourable directions of spins 
as the weakest scale. Our model treats the magnetic states in chiral magnets within a 
phenomenological quasi-classical continuum approximation. This micromagnetic approach 
is appropriate for chiral magnets where the helical modulations are large as compared to 
atomic distances, i.e., it applies to modulations with typical lengths of a few hundred up to 
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several thousand A observed in many chiral helix magnets 
specifically addresses intermediate length scales. 

We write the magnetic free energy density in the form 




. The model, hence, 



/ = Am" Y,^din,f + r^Aiymf + ^(m) + /o(m) 



(1) 



where the first and second term in rj and A {rj^A > 0), respectively, describe the mag- 
netic stiffness with n being a unit vector along the magnetization, m = mn). The 
third term //)(m) = DC{m) in the free energy describes the chiral interactions, where 
the Dzyaloshinsky-Moriya constant D determines its handedness and strength. The DM- 
interactions can be expressed in terms of so-called Lifshitz invariants, i.e., linear, antisym- 
metric gradient terms of the magnetization = [niidkmj — nijdkmi) where i,j,k are 
combinations of Cartesian coordinates x, y, z that are consistent with the symmetry class of 
the system under consideration. The fourth term foim), finally, includes nongradient terms 
of the free energy and may be expanded according to Landau theory in even powers of m. 
We neglect weaker contributions such as magnetic anisotropics, or dipolar interactions and 
address their importance in the discussion of our results. 

The magnetic stiffness terms in rj and A are motivated as follows. It is convenient to con- 
sider at first local moments, characterised by a fixed magnetisation modulus. For isotropic 
local moment ferromagnets (LMFM) in the ordered state it is well-known ^^l, that the 
magnetisation as described in a continuum limit of the Heisenberg model displays changes 
of orientation and modulus due to magnon-magnon scattering. This situation is described 
by the gradient energy given above for 77 = 1. Patashinskii and Pokrovskii IS] further 
showed that magnon-magnon scattering in LMFM also forces the longitudinal stiffness to 
vanish in the limit of infinite wavelength, i.e., rj = 1 is valid in this limit also (see also, 
Zwerger ^^). However, the SU(2) symmetry of the spin system enforces that the lon- 
gitudinal stiffness is the result of higher order transverse processes. Further mechanisms 
on top of magnon-magnon scattering, such as coupling to the particle-hole continuum in 
itinerant-electron systems, dipolar interactions and defects, therefore additionally reduce 
the longitudinal stiffness before it vanishes altogether in the limit of infinite wavelength. A 
minimal description of this additional longitudinal softening requires the second gradient 
term written in rj with 77 < 1, where t] expresses the ratio between longitudinal and trans- 
verse stiffness of the ordered magnetic state. For an extensive discussion of the technical 
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validity of the two gradient terms used here and their relationship to conventional spin fluc- 
tuation theories of itinerant-electron magnets, we refer to the supplementary information of 
this paper 



2C I published data of the longitudinal susceptibility in the local- 



We have also analysed 

moment ferromagnet EuS 21 1 and the itinerant-electron ferromagnets Ni [22] and MnSi [2^ . 
Our analysis establishes values of rj for EuS, Ni and MnSi of 0.925, 0.65 and 0.4, respec- 
tively, i.e, important reductions of the longitudinal stiffness exist in all real materials. The 
experimental data for these archetypical ferromagnets show that the regime 77 < 1 in our 
model is valid and relevant for many different types of magnetic systems. In particular, the 
very prominent longitudinal softening in itinerant-electron magnets Ni and MnSi has not 
yet been appreciated before. In the presence of chiral interactions the reductions of longi- 
tudinal stiffness become of tantamount importance, because they facilitate the formation of 
spontaneous skyrmion ground states. 

For our model we first show that certain individual skyrmions represent localised re- 
gions of reduced ground state energy on a background of vanishing magnetization. General 
consideration of the stability of particle-like solutions as function of dimension against com- 
pression/stretching |12] strongly suggest that two-dimensional solutions are the most stable. 
We therefore begin by considering an individual skyrmion of cylindrical symmetry, as shown 
in Fig. 1 b and c to explain the basic physical ideas. As explained in further detail below, the 
stability of such a cylindrical structure establishes the possibility of spontaneous skyrmion 
lattices in all noncentrosymmetric crystallographic classes that allow chiral interactions Q| . 

Expanding the contribution fo{m) according to Landau theory for small m as /o = a{T — 
Tc)m? + hrn^ -|- . . .. leads in the absence of chiral interactions (/d = 0) to the conventional 
Curie temperature Tc of centrosymmetric systems. Thus, when the temperature drops 
below Tc the energy density of a ferromagnetically spin-aligned state is lowest. Chiral 
interactions {fo < 0) favouring rotations of the moments with respect to each other reduce 
the energy density further. This stabilization of twisted magnetic states occurs only by the 
competition of the DM-interactions with the exchange, expressed by the form of the gradient 
energy. Therefore, the exact form of /o(^) is not decisive for the stability of skyrmions. In 
particular a form of /o(^) including higher order terms with 6 < and a stabilizing term 
for a first-order transition will change only quantitative relations and the details of the 
phase-diagram. Possibly, it may lead to the discontinuous nucleation of skyrmions. Figure 1 
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illustrates the twisted magnetic structures with rotating magnetization and the associated 
energy densities. For rotations along a single given direction, shown in Fig. 1 a, the reduction 
in energy density, shown in Fig 1 d, is uniform. The transition to the ferromagnetic state 
is consequently pre-empted by a transition to a one-dimensionally modulated state at a 
temperature 

As shown in Fig. 1 b and d rotations of the magnetisation in two dimensions reduce the 
energy density even further near its nucleation point at the center of the cylindrical structure 
p = 0. Considering only the energy reduction at p = we obtain as an upper limit for the 
temperature. To, for the onset of this so-called double-twisting 

The characteristic temperatures T/j and depend only on a few physically meaningful, 
material specific parameters that may be determined from experiment. The ratio A/D 
measures the pitch of the chiral modulations in the ordered helical state, while D /a measures 
the chirality of magnetic fluctuations in the paramagnetic state. Hence, D/a is related to a 
chiral component of the susceptibility. The parameters a and b are the initial susceptibility 
and mode coupling parameter, respectively. This permits estimates of Th and To from 
experiments. An example is given below. 

As shown in Fig. 1 d the energy reduction due to a double-twist is no longer uniform. 
In fact, for fixed magnetisation amplitude the energy reduction turns into an excess of the 
energy density at larger distances p from the center, which outweighs the initial reduction. 
This destabilises the double-twist structure altogether, which is the reason why magnetic 
skyrmions prior to the work reported here were believed to be unstable. Skyrmion textures 
may only form spontaneously when no additional energy is required for double-twist rota- 
tions at large p. The latter may be achieved by permitting the magnitude of the moment 
to be soft. As shown in Fig. Ic and Id skyrmions under these conditions have the lowest 
energy density, when the moment decreases with increasing distance from the centre. 

The results of a comprehensive numerical analysis of the free energy for skyrmions taking 
into account the variation of the energy density with p for differing rj are summarised in Fig. 2 
and 3 (see Eqn. ((T)) and Methods section). This analysis establishes rigorously that stable 
cylindrical skyrmions form spontaneously in a finite temperature interval Tt < T < Ts oi the 



phase diagram Fig. 2, where the hmits Th and are also marked. The magnetic structure of 
the skyrmions is characterised by a magnetisation vector m = (m, 6', ip) that rotates outward 
from the axis in all directions, while the modulus m(p) shrinks and decreases towards the 
boundary of the tube at radius R (Fig. Ic). The amplitude of the local magnetisation 
m(p) vanishes continuously when approaching T5 from below, as shown in Fig. 3. Thus, 
the skyrmion nucleation and the formation of the skyrmion texture coincide. The increase 
of the magnetization with decreasing temperature drives the transition from the skyrmion 
phase to the uniform helix phase at a lower temperature T^. The transition into a one- 
dimensionally modulated state by transforming the structure Fig. 1 b or c into the helix 
Fig 1 a cannot be performed by a smooth deformation of the skyrmions, because these two 
textures have different topology. For this reason, a first-order transition is generally expected 
at a temperature Tj. In the phase diagram, the transition line Tj has been determined from 
the approximate condition that the average energy of the circular skyrmion cell is equal to 
that of the helix. The skyrmion phase vanishes as T5 — Tt in the limit 77 — 1. This feature 
of the phase diagram (Fig. 2) concurs with an analytical result by Wright and Mermin for 
magnets with the properties of fixed moment classical Heisenberg models (Eqn. with 
7] = 1) j^. The result was believed to prove the absence of skyrmion ground states in chiral 
ferromagnets. However, it applies only to a particular form of the model where = 1, while 
the regime r/ < 1 had not been considered systematically. 

Having established that cylindrical skyrmions with soft amplitude can be localised regions 
of reduced energy we next address the proliferation of these skyrmions to form textures. In 
a chiral magnet, the formation of an extended skyrmion condensate is subject to: (i) weak 
residual magnetic anisotropics and dipolar interactions in the crystalline environment, and 
(ii) the magnetisation in the regime joining the skyrmion tubes. Under the most favourable 
conditions, the skyrmion condensate may form a regular lattice in analogy to the vortex 
lattice of type 2 superconductors. In fact, there exists a fundamental analogy between mag- 
netic skyrmion lattices and Abrikosov lattices. The formation of stable spatially modulated 
states in superconductors and in magnets can be justified by the concept of negative domain 
wall energies We illustrate the stability of a skyrmion condensate for this case 

by the full solution for an extended state in an isotropic system in two dimensions. We 
note that the solution as shown in Fig. 4 applies also to three-dimensional systems, when 
the solutions are artificially constrained to remain homogeneous in the third spatial direc- 
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tion (perpendicular to the plane shown in Fig. 4). The structure of this solution shows in 
particular that the core region of the skyrmions is unaffected, as stated above, while the 
magnetization between the skyrmions mediating their interactions does not even need to be 
strictly zero for the skyrmion lattice to become the ground state. 

The spontaneous formation of skyrmions based on cylinder-symmetric solutions applies 
to all symmetry classes including higher symmetries, say the cubic and the isotropic case 
without inversion center, because the core structure of the skyrmions is determined by the 
two largest energy scales, namely the exchange energy and chiral interactions, while the 
texture is determined by weaker interactions. Technically the general validity may be seen 
by noting that the Euler equations for the modulus of the magnetisation m and the polar 
angle 6 which control the energetic stability of the skyrmions, given by Eqn. (7) and (8) 
in the Methods section, respectively, are identical for all relevant symmetry classes (see 
Refs.lllQ for a discussion). In contrast, the behaviour of the azimuthal angle ip, which 
differs for the different symmetry classes, only affects the manner of rotation of the magnetic 
moments in the radial direction. General considerations of the stability of localized textures 
br chiral magnets suggest that droplet-like skyrmions are not stable in three dimensions 



12j | . Therefore, we believe that the two-dimensionally localized skyrmion tubes discussed so 



far will form extended ground-state condensates in three-dimensional systems. However, in 
a three-dimensional environment bending and terminating skyrmion tubes imposes a cost 
in energy that determines the nature of the condensate. The additional cost in energies 
may be readily reduced by closing the tubes amongst themselves or ending them at surface 
and defects. In isotropic three-dimensional materials the skyrmion condensate will therefore 
have a liquid-like or amorphous appearance driven by the frustration between preferred 
orientations and their elastic interactions in the condensed structure. 

Experimental systems where skyrmion ground states may exist will have to display soft- 
ened longitudinal stiffness together with a hirachy of energy scales in which spin-aligning 
exchange interactions are dominant followed by subleading chiral interactions causing spin 
rotations. Further, magnetic anisotropics must be low because they prevent the continuous 



39|. The el- 



rotation of the magnetization and, thus, multiply-modulated textures jl^ . IgI 
ementary mechanism, shown here for a ferromagnet and illustrated in Fig. 1, then applies 
to any chiral system with a tendency to form twisted structures. In particular, our model 
also applies to antiferromagnetic crystals with chiral symmetry by replacing the magnetiza- 
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tion with a staggered antiferromagnetic vector order parameter 38|. Owing to the broken 
inversion symmetry at the surfaces of magnetic materials l^d | thin films and interfaces 
are the largest class of condensed matter systems where skyrmion condensates may form 
spontaneously. Appropriate metallic films of this kind have not been investigated in detail 
in the regime of interest. A less frequent, but equally promising class of materials are bulk 
materials with chiral interactions due to a lack of inversion symmetry. 

The most powerful technique to test our prediction experimentally is polarised magnetic 
neutron scattering. For a periodic skyrmion lattice neutron scattering will show long-range 
order akin to the vortex lattice in superconductors, while the scatteririg intensitv of an 
amorphous appearance will be akin to partial order in liquid crystals 40|. In the bulk 
properties subtle evidence for the predicted intermediate phase is expected, when going 
from paramagnetism to helical magnetic order. Real-space imaging of magnetic textures by 
magnetic force microscopy, spin sensitive tunneling, magnetic x-ray microscopy, and Lorentz 
transmission electron microscopy are finally expected to provide evidence of skyrmion lattices 
akin vortex lattices in the Shubnikov phase of superconductors. The helical structure in 
the chiral magnet (Fe,Co)Si were recently imaged directly for the first time j3| because 
of the exceptionally long helix period in this material, but the technique may introduce 
uncontrolled strains that interfere with the intrinsic properties. It is to be expected that 
the resolution of improved imaging techniques will soon be sufficient to take pictures of the 
skyrmions predicted here. 

Our study was inspired by recent high pressure experiments in the cubic itinerant-electron 
helimagnet MnSi. The low temperature resistivity in MnSi changes abruptly from a T^-Fermi 
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liquid temperature dependence to a T^/^ dependence above a pressure Pc = 14.6 kbar 
The resistivity is unchanged up to at least 3pc, suggesting the formation of a stable phase. 
However, as emphasized in Ref. j33] a T'^/^ resistivity is normally characteristic of spin glasses 
and amorphous ferromagnets where the conduction electrons undergo a diffusive motion, but 
not for the ultra-pure MnSi single- crystals investigated, which by all accounts should show 
the behaviour of a Fermi liquid. The spontaneous formation of an amorphous magnetic state 
in ultrapure samples would provide a simple explanation that resolves this contradiction of 
high purity and diffusive charge carrier motion in MnSi. It also explains the stability of the 
resistivity as a new phase. Magnetic neutron scattering in MnSi further revealed, that large 
ordered moments are present in the nominally paramagnetic phase. The ordered moments 
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are organised such that they lead to scattering intensity on the surface of a sphere with 
unexplained broad maxima for (110). If we suppose that the easy magnetic axis, (111), does 
not change as function of pressure, while the magnetism changes from a helical modulation 
along (111) to cylindrical skyrmion tubes with their axes preferentially along (111) we expect 
intensity in the great circles perpendicular to (111). These great circles interset along the 
(110) directions, thus explaining the broad maxima in this location. Finally, as shown 
above we predict a spontaneous skyrmion phase prior to the onset of helical order. From 
the values of Tq = 29.5 K, the wavelength of the helix go = 0.039 and the effective 
Curie- Weiss moment pejf = 2.2jiB we estimate To — Tc^ 0.9K [2^. A hump in the specific 
heat above Tc and and an associated sphere in neutron scattering intensity as observed in 
ambient pressure experiments in MnSi are consistent with such a phase over a temperature 
interval of ~ 1 K |20l . l3l|. We note that recent theories on the partial magnetic order in 



MnSi by Tewari et al. and Binz et al. 



27| also suggest inhomogeneous twisted states 



with modulated amplitude. These theories rely on unsubstantiated higher-order interaction 
terms and, therefore, share the problems with the corresponding non-linear field theories 
without identifying the basic mechanism of stable skyrmions in such textures. The evidence 
for a skyrmion ground state in MnSi discussed here resolves all of the existing contradictions, 
but a definitive proof will require a host of novel techniques listed above. 

We finally note, that the argument for the existence of skyrmion textures in ferromagnetic 
metals is closely related to the blue phases of chiral nematic crystals 0,0]. The standard 
continuum model of chiral nematics admits localized structures as vortex tubes [sS], where 
the energy disadvantage by double- twisting a configuration has been recognised early on . 
However, in the blue phases topological defects in the form of disclination lines, which do 
not exist in magnetic systems, allow to overcome the disadvantage in energy and stabilise 
extended textures. In chiral magnets skyrmion textures appear instead as condensate of 
strictly localized skyrmions with lines or sheets of vanishing magnetic order between them 
(Figure 4). Fluctuations and the competition betwen different frustrated couplings entail 
the possibility of internal transitions between various textures with liquid, glassy, or lattice- 
like organization of skyrmions. Thus, the predicted magnetic skyrmion textures promise 
rich phase diagrams, similar to those of liquid crystals and Abrikosov lattices of flux-lines 
in type-II superconductors Isd j. 
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Methods 



Cubic noncentrosymmetric magnets (e.g., MnSi and other B20 compounds |15|) ) belong 
to the crystallographic class T. For this class the Dzyaloshinsky-Moriya energy is 



fD = D (4;) + C^y} + 4^)) = Dm ■ rotm, 



zy 



(4) 



where D is the Dzyaloshinsky constant [15']. The helix solutions for these systems are 
described by a linear wave with a wave number an = D I {2 A) and constant magnetisation 
am plit ude rri^ = 2a{Th — T)/b j9|, |15|- For all uniaxial systems the invariants are listed 
in |l6|. In particular, for uniaxal magnets with Cnv symmetry {n = 3, 4, 6) one has = 
D{Ct^ - Ci^). The n-fold axis is directed along z. Examples are Tb3Al2 or Dy3Al2 with 
tetragonal noncentrosymmetric structure [14] or the superconducting antiferromagnet 
CePtaSi, space group P4mm j4l]]. The structure of the cylindrically symmetric skyrmions 
is given by solutions of the Euler equations for the free energy Eqn. (^. The magnetization 
vector m = (m, 6, ip) is aligned along the cylinder or 2;-axis in the skyrmion center, p = 0. 
In cubic and uniaxial systems solutions for a skyrmion homogeneous along the 2;-direction 
are given by m(p), 0{p), ip = sign0 + 0o, where sign = + or — and (po is a constant, both 
determined by the specific symmetry class The dependence m(p) and 6{p) describes 
the variation of the amplitude of the magnetic moments and their angle as function of the 
distance p to the cylinder axis in cylindrical coordinates, respectively. After substitution of 
the solutions for il'{4') and an integration with respect to z and the energy of the skyrmion 
is reduced to the following functional jFg = 2nL /s(^, 6)pdp with energy density 



/rf^y ^ sin^^ 



de 

dp 



+ 



sin 9 cos 9 



P 



(5) 
(6) 



\dp J p^ 

where L is the length of the skyrmion along ^-direction, and the spatial coordinate p is 
measured in q^^ units. The Euler equations for the functional © 



d'^m 
1^ 



1 / dm\ 
p\dp J _ 



Agi9,9,)m-q^^^ = 



(7) 



m 



£9_ 

dp^ 



1 (d9\ sin cos 2sin^9 



p \dp J 



P^ 



P 



+ 2 



d9 ^ 
dp 



d'm\ ^ 
dp J 
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with the boundary conditions ^(0) = 0, m(0) = mo, m{R) = describe the structure of 
the isolated skyrmion and yield the energy as a function of mo and R. The equilibrium 
skyrmion profiles are derived by minimization of ^^^(mo) with respect to mo and R; the 
radially integrated energy densities of the vortices per lengths (Figure 1 d) are given by 
w{p) = fs{m,9)pdp. Full solutions for extended two-dimensional models have 

been constructed numerically using a finite-difference scheme on rectangular lattices with 
dynamically optimized lattice parameters and periodic boundary conditions. The search for 
ground-states and checks on their stability were done using a continuous variable Monte 
Carlo simulated annealing method and mesh- refinement. 

The unit length used throughout this paper is l/qo, which is the inverse of the twisting 
length in chiral magnets with Dzyaloshinsky-Moriya couplings. Energy and temperature are 
given in equivalent units of D'^/{2Aa) as in Eqs. © and ©■ 
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FIG. 1: Chiral modulated structures for noncentrosymmetric ferromagnets (a - c). a, one- 
dimensional modulated structure propagating along a direction x and with a constant modulus of 
magnetic moments, b, Circular cross section through a skyrmionic excitation with fixed modulus. 
It can be visualized by rolling one half-period of a helix modulation as in a into a cylindrical tube 
with propagation axes now pointing radially out in all directions p from the center, c, A circular 
skyrmion cell. The magnetization has maximum modulus at the skyrmion axis and rotates along 
the radial directions {ip = (p). The modulus m gradually decreases with increasing distance from the 
center, and equals zero at the boundaries of the cell. The chirality (handedness) of the cylindrical 
skyrmions shown here may be seen as follows: when the fingers on a hand follow the moments rolling 
outward the thumb defines a winding sense (chirality) around the center (note that space inversion is 
not satisfied with respect to points in the center as this winding direction would have to be reversed 
upon space inversion). The magnetic structures a, b, c can exist in noncentrosymmetric crystals 
of class Cnv The chiral rnagnets Tb3Al2 or Dy3Al2 with tetragonal noncentrosymmetric structure 
C|j, belong to this class [ij]. For this crystallographic class the magnetization rotates around an 
axis, that is perpendicular to the propagation directions in the one-dimensionally modulated state 
and tangentially for the cylindrical skyrmions. In noncentrosymmetric crystals from the classes 
D„ and in cubic or isotropic systems, the magnetization in the helix and the skyrmions rotates 
in the planes perpendicular to the propagation directions, i.e., about the propgagation directions. 
Such structures can exist in the noncentrosymmetric ferromagnets with B20-structure, in particular 
in MnSi 1 15]. For other crystallographic symmetries more involved structures can occur. All of 
these magnetic structures in three-dimensional crystals can be described either as one-dimensional 
modulated state, which is a helix, or as cylindrical skyrmions. Phenomenological theory allows 
to analyse the stability of these structures for different crystallographic classes within a common 
mathematical framework, see Methods section and Ref. 'l(y\. d, comparison for the local energy 
density per unit volume for the helical modulation, a, and fixed and varying magnetic moment 
m, b and c, respectively, in circular skyrmion cells (radius R) along a radial direction p. The 
corresponding helical state has a constant energy density which is always higher than the energy 
density in the center of the skyrmions. For fixed-moment skyrmions this energy gain is offset by a 
large excess energy near the boundary of the cell at R. 
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FIG. 2: Phase diagram of a chiral ferromagnet in terms of temperature versus longitudinal stiffness- 
parameter {rj,T). The diagram is based on the stability of the cylindrical skyrmions shown in 
Fig. 1 c. Therefore, the diagram applies for two- and three-dimensions. Above the ferromagnetic 
Curie temperature Tc two characteristic temperatures, defined in Eqn. 0, rule the phase diagram: 
Td for the instability against double twisting and Th for the formation of the helical phase. The 
skyrmion phase is thermodynamically stable between the line Ts{ri) for the continuous transition 
into the paramagnetic phase and the line Tt{r]) for the first-order transition into the helix. This 
line is independent on the chosen parameter b. The helical phase exists as a metastable phase 
between and Tt{rj). The estimate for the first-order line Tt{r]) has been calculated for models 
with h = 0.05. 
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FIG. 3: Details of skyrmions solutions for rj = 0.4 and b = 0.05 at various temperatures. (See, 
Eqs. (O,® in Methods section.) a, Variation of the modulus of the magnetisation as function 
of radius at different temperatures, b, Variation of the angle of rotation, 6, as function of radius 
at different temperatures. At low temperatures the magnetization rotates from angle zero in the 
center to angles 6{R) < vr at the cell boundaries. The angles 6{R) depend on temperature. The 
radius is given in units of the twisting lengths l/qo- c. Evolution of the radius of the skyrmion 
cell up to the transition into the paramagnetic state. In the transition region both the radius of 
the skyrmion cell R and the angles 0(R) grow unlimitedly by approaching the critical temperature 
Ts- d, Corresponding evolution for the amplitude magnetization in the center m(0). The modulus 
gradually decreases to zero as the temperature approaches the critical value Ts while the radius 
of the skyrmions diverges. By this process the single skyrmions forming a dense stable condensate 
below Ts can transform into isolated skyrmions. Note that the characteristic twisting length 
in radial direction, which is given by the slope 0{p) in panel b, does not change with changing 
temperature. The core region of the skyrmions with sizeable magnetization has essentially a fixed 
diameter, while only a tail with diverging radius and rapidly vanishing local magnetization evolves 
near Ts- This diameter is given by the ratio D/A, which also fixes the characteristic size of chiral 
fluctuations in the paramagnetic phase above Ts- 
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FIG. 4: Structure of a two-dimensional skyrmion lattice, derived as minimum energy solution for 
the model Eqn. with Dzyaloshinskii-Moriya interactions, Eqn. Q. The ground-state has four- 
fold symmetric lattice structure, and net magnetization direction of the skyrmion cores is staggered. 
Thus, the periodicity of this square lattice is doubled as a 2x2 antiferromagnetic structure, a, 
overview showing the modulus m and the corresponding magnetization vectors, b, Projection 
of the magnetization vectors into the base-plane. Units shown in b are the twisting length I/qq. 
Red/blue signals local magnetization direction with positive/negative components arbitrarily taken 
perpendicular to the two-dimensional plane. For the core-regions of the skyrmions, lines of constant 
modulus, m =const, are shown as contour lines in the base-plane. This result applies to thin 
magnetic films with broken inversion symmetry made from isotropic or cubic (crystallographic 
class T) metallic materials. The magnetization direction rotates around the radial directions going 
outward from the skyrmion center. A three-dimensional analogue of this dense texture composed 



of an amorphous arrangement of cylindrical s 



in the cubic B20 metals FeGe and MnSi 



3C, 



cyrmion strings is consistent with recent experiments 
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as discussed in the text. 
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